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The papain-like protease (PL°’) of Middle-East respiratory syndrome coronavirus (MERS-CoV) has 
proteolytic, deubiquitinating, and delSGylating activities. The latter two are involved in the 
suppression of the antiviral innate immune response of the host cell. To contribute to an 
understanding of this process, we present here the X-ray crystal structure of a complex between 
MERS-CoV PL®” and human ubiquitin (Ub) that is devoid of any covalent linkage between the two 
proteins. Five regions of the PL” bind to two areas of the Ub. The C-terminal five residues of Ub, 
RLRGG, are similar to the P5—P1 residues of the polyprotein substrates of the PL” and are 
responsible for the major part of the interaction between the two macromolecules. Through site- 
directed mutagenesis, we demonstrate that conserved Asp165 and non-conserved Asp164 are 
important for the catalytic activities of MERS-CoV PL”. The enzyme appears not to be optimized 
for catalytic efficiency; thus, replacement of Phe269 by Tyr leads to increased peptidolytic and 
deubiquitinating activities. Ubiquitin binding by MERS-CoV PL” involves remarkable differences 
compared to the corresponding complex with SARS-CoV PL””. The structure and the mutational 
study help understand common and unique features of the deubiquitinating activity of MERS-CoV 
pro 
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INTRODUCTION piratory syndrome coronavirus (SARS-CoV; Drosten et 


al., 2003; Ksiazek et al., 2003; Kuiken et al., 2003; Peiris 


To date, six coronaviruses infecting humans have been 
characterized. Infections with human coronaviruses 
(HCoVs) 229E (Hamre & Procknow, 1966), OC43 
(McIntosh et al., 1967), NL63 (van der Hoek et al., 
2004), and HKU1 (Woo et al., 2005) cause relatively 
mild symptoms in most cases, whereas severe acute res- 
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et al., 2003) and Middle-East respiratory syndrome 
coronavirus (MERS-CoV; Zaki et al., 2012) are connec- 
ted with severe respiratory-tract infection and, in particu- 
lar in case of MERS-CoV, acute renal failure (Eckerle et 
al., 2013), leading to high case-fatality rates of ~10 and 
~35%, respectively. In spite of 13 years of research on 
SARS-CoV (Hilgenfeld & Peiris, 2013), no approved 
drugs or vaccines are available for the treatment or pre- 
vention of coronavirus infection (Wang et al., 2016). 
This is mainly due to the fact that although these emer- 
ging viruses have devastating effects on those infected, 
the absolute numbers of cases (~8000 for SARS, 1733 so 
far for MERS; (http://www.who.int)) imply that the de- 
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velopment of specific antivirals is very likely not com- 
mercially viable. On the other hand, the global risk posed 
by MERS-CoV must not be underestimated. Since its 
discovery in September 2012, the number of MERS 
cases reported has been rising steadily, with some inter- 
mittent peaks connected to hospital outbreaks in Saudi 
Arabia (Assiri et al., 2013). Man-to-man transmission of 
MERS-CoV has also been impressively demonstrated by 
the recent outbreak of MERS in South Korea, which was 
traced back to a single traveller from the Arab peninsula 
(Butler, 2015). Therefore, it is imperative that academic 
laboratories help increase the preparedness against a pos- 
sible MERS-CoV pandemic by characterizing antiviral 
drug targets and by identifying lead compounds interfer- 
ing with them. 

In order to successfully infect humans, a virus has to 
meet at least two conditions: 1), it should maintain a suf- 
ficiently correct replication of its genetic material; 2), it 
should inhibit the host antiviral response. The papain- 
like protease (PL””’) of MERS-CoV (or SARS-CoV) is 
involved in both of these tasks (Yang et al., 2013; Bar- 
retto et al, 2005). The PL” is a domain located in the 
middle part of the largest non-structural protein, Nsp3, of 
MERS-CoV (or SARS-CoV). It is responsible for releas- 
ing Nspl, Nsp2, and Nsp3 from the polyproteins la 
(ppla) and lab (pplab), an essential step of replication 
(Harcourt et al., 2004). Like its SARS-CoV counterpart, 
the MERS-CoV PL”” also has deubiquitinating (DUB) 
and deISGylating activities in vivo as well as in vitro 
(Yang et al., 2013; Mielech et al., 2014; Lei et al., 2014; 
Baez-Santos et al., 2014b). K48- and K63-polyubiquitin 
poly (Ub) and ISG15 (interferon-stimulated gene 15) - 
conjugated targets are usually involved in host innate im- 
mune regulation (Liu et al., 2005; Maringer & Fernan- 
dez-Sesma, 2014). The PL”” has the ability to digest 
K48- and K63-linked polyUb chains and to remove 
ISG15 from ISG15-linked proteins (Baez-Santos et al., 
2014b), thereby interrupting the signalling pathways 
leading to the innate immune response. Thus, the PL”” 
can block the activation of IFN regulatory factor 3 
(IRF3) (Yang et al., 2013) and subsequently the produc- 
tion of interferon B (IFNB) (Mielech et al., 2014). Inter- 
estingly, MERS-CoV PL”” shows a similar cleavage rate 
for K48- and K63-linked polyUb chains, while the 
SARS-CoV enzyme prefers K48- over K63-linked 
chains (Baez-Santos et al., 2014b). The former enzyme 
degrades a polyUb chain by removing mono-Ubs, where- 
as the latter cleaves di-Ub units off the polyUb chain 
(Bekes et al., 2015). 

We have reported the first crystal structure of the 
MERS-CoV PL” (Lei et al., 2014). Later, two other 
groups also described the structure of this enzyme 
(Bailey-Elkin et al., 2014; Lee et al., 2015). The struc- 
ture of PL” can be divided into two parts: a ubiquitin- 


like (Ubl) domain and a catalytic domain with thumb, 
palm, and fingers subdomains. The overall fold of 
MERS-CoV PL?” is not only similar to that of SARS- 
CoV PL”, but also to that of several human ubiquitin- 
specific proteases (USPs) (Hu et al., 2005). In 2014, the 
X-ray structure of the complex of SARS-CoV PL” with 
ubiquitin has been reported (Chou et al., 2014; Ratia et 
al., 2014). Several key residues (such as Glu168 or 
Tyr265) of SARS-CoV PL”” that are important for ubi- 
quitin recognition (Chou et al., 2014; Ratia et al., 2014), 
are not conserved in MERS-CoV PL”. Bailey-Elkin et 
al. (2014) described the structure of an artificially linked, 
covalent complex between ubiquitin and MERS-CoV 
PL””. Here, we present the crystal structure of a non-co- 
valent complex between the two proteins and a mutation- 
al study of the interactions involved. For these studies, 
we used the Cys111Ser active-site variant of MERS- 
CoV PL”. 


MATERIALS AND METHODS 


Recombinant production of MERS-CoV PL” 
and its variants 

The PL’ of MERS-CoV (strain 2c EMC/2012; Gen- 
Bank no. AFV09327.1) contains 320 residues, from 
Gln1482 to Asp1801 of ppla/lab. For simplification, GIn1482 
was renumbered into Gln! here. The DNA plasmid cod- 
ing for MERS-CoV PL” was produced earlier (Lei et 
al., 2014). 

The MERS-CoV PL” C111S, D164E, D164A, D165E, 
D165A, and F269Y variants were produced using the 
same strategy that we described before (Lei et al., 2014). 
All primers for these variants are listed in Supplemental 
Table S1. All DNA plasmids coding for the altered PL”” 
were verified by sequencing. 

Genes coding for wild-type (WT) MERS-CoV PL’” 
and for its variants were expressed and the correspond- 
ing proteins were purified according to our previous de- 
scription (Lei et al., 2014). 


Recombinant production of SARS-CoV PL” 

The PL’ of SARS-CoV (Strain TOR2; GenBank no. 
AY274119.3) comprises 319 amino-acid residues, cor- 
responding to Glul1541 to Tyr1859 of ppla/lab. A gene 
coding for the SARS-CoV PL”™°was amplified by PCR 
with the following two primers 5'"-CTAGCTAGCGAGGTTA 
AGACTATAAAAGTGTTC-3' (forward) and 5'- 
CCGCTCGAGTTAATACGACACAGGCTTGATGGTTG 
TAG-3' (reverse). The PCR product was digested by Nhe 
I and Xho I, then was ligated into the pET-28a plasmid 
(Novagen). The recombinant plasmid DNA was verified 
by sequencing. Expression of the gene construct coding 
for SARS-CoV PL” and purification of the protein were 
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performed according to the procedure described for 
MERS-CoV PL” (Lei et al., 2014). 


Crystallization of MERS-CoV PL”” with ubiquitin 
Purified MERS-CoV PL’”® (WT) and a variant that had 
the active-site Cys111 replaced by Ser (PL”” (C111S)) 
were both concentrated to ~24 mg/mL in 20 mmol/L 
Tris-HCl, 150 mmol/L NaCl, pH 8.8, 10 mmol/L B-mer- 
captoethanol (BME). Human ubiquitin (BostonBiochem) 
was dissolved to 6 mg/mL in 20 mmol/L Tris-HCl, 150 
mmol/L NaCl, pH 8.8. 500 wL PL” (WT) or PL” (C111S) 
were mixed with 500 wL ubiquitin (~1:1 molar ratio) at 4 
°C overnight. The complex of PL” (WT)-Ub or PL” 
(C111S)-Ub was purified by gel filtration (HiLoad™ 
16/60 S200 column, GE Healthcare) in 20 mmol/L Tris- 
HCl, 150 mmol/L NaCl, pH 8.8 the next day. The final 
concentration of PL”” (WT)-Ub or PL?” (C111S)-Ub 
was ~12 mg/mL. The two complexes were crystallized 
using the sitting-drop method and a Phoenix crystalliza- 
tion robot (Art Robbins) at 18 °C. 0.25 uwL of protein and 
0.25 uL of reservoir were mixed and equilibrated against 
75 pL reservoir. Screening kits Index’, SaltRx™, PEG 
Rx 1 & 2, PEG/Ion © 1 & 2 (Hampton Research), and 
Structure Screen 1 MD1-01, Structure Screen 2 MD1-02 
(Molecular Dimensions) were used. Crystals of PL”” 
(C111S)-Ub were observed under condition No. 9 of 
MD1-01, whereas no crystal of PL” (WT)-Ub was ob- 
tained. Optimized crystals of PL” (C111S)-Ub were 
subsequently obtained within two days under the condi- 
tion: 22% w/v PEG 4000, 15% v/v 2-propanol, 0.1 mol/L 
tri-sodium citrate pH 4.8, and 10% glycerol. 2 uL of pro- 
tein and 2.5 wL of reservoir were mixed to equilibrate 
against 500 pL reservoir. 

Crystals were placed in a nitrogen-gas stream (100 K). 
A 3.16-A dataset was collected at wavelength 0.91841 A 
at beamline 14.2 of BESSY, Berlin (Mueller et al., 
2012). The program XDS (Kabsch, 2010) was used to 
process the diffraction data. The space group was found 
to be P6;, with unit-cell parameters a = b =138.14 A, c= 
57.59 A, y = 120°. Diffraction data statistics are shown 
in Table 1. 


Phase determination, model building and 
refinement 

The structure of the MERS-CoV PL” (C111S)-Ub com- 
plex was solved by molecular replacement using 
MOLREP (Vagin & Teplyakov, 2010). The program se- 
lected the MERS-CoV PL?” (Protein Data Bank (PDB) 
entry 4P16, Lei et al., 2014) as the first search model. 
Human ubiquitin (PDB entry: 1UBQ, Vijay-Kumar et 
al., 1987) was used as the second search model. The 
model of the complex was inspected and rebuilt using 
Coot (Emsley et al., 2010), and refined using phenix.re- 
fine (Headd et al., 2012). The refinement statistics are 
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Table 1. Data collection and refinement statistics 


MERS-CoV PL*® (C111S)-Ub 


Data collection statistics 
Space group 
Unit-cell dimensions (A, °) 


Wavelength (A) 

V,, (A°/Da) 

Solvent content (%) 
Resolution range (A) 
Number of unique reflections 


P6, 

a=b=138.14, c=57.59 
y=120 

0.91841 

3.59 

65.8 

45.22-3.16 (3.33-3.16) 
10955 (1577) 


Roe 0.050 (0.413) 
CC(1/2) 0.997 (0.708) 
Completeness (%) 99.9 (100) 
Mean I/o (1) 12.8 (2.0) 
Multiplicity 7.3 (7.5) 
Refinement statistics 

ee COP pa 

Riree (%)° 25.5 

No. of atoms 

Protein 3054 

lon 1 

Ligand 24 
Clashscore* 17 
r.m.s.deviation in bond lengths (A) 0.012 
r.m.s.deviation in bond angles (°) 1.945 
Average B-factor for all atoms (A?) 108 
Average B-factor for PL”® (A’) 95 
Average B-factor for Ub (A’) 145 
Ramachandran plot 

Residues in favored regions (%) 95.30 
Residues in allowed regions (%) 4.44 
Residues in outlier regions (%) 0.26 


Notes: ' Rpim (Weiss & Hilgenfeld, 1997). * Roryst = Dnxi 


| Fo(hkl)-F (hk!) | dig F (kl). Rie WaS Calculated for a test 
set of reflections (4.9%) omitted from the refinement. 

* Cl-ash-score is defined as the number of clashes calculated 
for the model per 1000 atoms (including hydrogens) of the 


model (Chen et al., 2010). 


shown in Table 1. The final model coordinates and struc- 
ture factors have been deposited in the PDB database 
with the code 4WUR. Figures (except for the supple- 
mental figure) have been prepared using Pymol (Schr6- 
dinger; http://www.pymolL.org/). 
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Kinetic assays of purified PL’’s 

All enzymatic assays were performed using a 96-well 
microtiter plate and the reaction buffer 20 mmol/L Tris- 
HCl, 150 mmol/L NaCl, pH 7.9, 2 mmol/L dithiothreitol 
(DTT). The fluorogenic substrates Cbz-Arg-Leu-Arg- 
Gly-Gly-7-amino-4-methylcoumarin (Z-RLRGG-AMC) 
(Bachem), Z-LRGG-AMC (BostonBiochem), and ubi- 
quitin- AMC (Ub-AMC) (BostonBiochem) were used. 
The fluorescence of free AMC with different concentra- 
tions (5 nmol/L—2.5 umol/L) in reaction buffer was 
measured to generate a calibration curve, in order to con- 
vert the change of fluorescence intensity per unit of time, 
A(AFU) /s, into the amount of hydrolyzed substrate in 
umol/L/s. 

The enzymatic cleavage reactions were run with an 
F1x800 fluorescence spectrophotometer (BioTek), to 
measure the increased fluorescence signal (A,,: 360 nm; 
hem: 460 nm) resulting from AMC release. Reactions 
were initiated by addition of the proteases to the reaction 
system. The peptide-hydrolysis kinetic assays were per- 
formed with the following conditions: 1 umol/L MERS- 
CoV PL” variant (D164E, D164A, D165E, F269Y), or 
10 pmol/L D165A, or 0.1 pmol/L SARS-CoV PL”, with 
different concentrations (10, 20, 40, 80, 100 umol/L) of 
Z-RLRGG-AMC or Z-LRGG-AMC in a final volume of 
100 pL at 25 °C. The kinetic curves for the proteases and 
their variants with the substrates Z-RLRGG-AMC or Z- 
LRGG-AMC were linear and the initial velocities also 
increased linearly with substrate concentration. No satur- 
ation was observed. Therefore, the data were fitted to the 
equation v/[E]}.o., = Kay)LS], where k,,,, approximates k_,/Kyy, 
as described previously (Barretto et al., 2005; Wojdyla et 
al., 2010). 

The deubiquitinating kinetic assays were performed 
under the following conditions: 0.1 wmol/L MERS-CoV 
PL’” wild-type or its variants D164E, D164A, D165E, 
F269Y, or 0.5 umol/L D165A, or 0.025 umol/L SARS- 
CoV PL” were incubated with increasing concentra- 
tions (1, 2, 4, 8 umol/L) of Ub-AMC in a final volume of 
50 uL, at 25 °C. Although PL” actually cleaves the iso- 
peptide bond between the carboxyl] group of the C-ter- 
minal Gly in Ub and the s-amino group of Lys in ubi- 
quitinated targets in vivo, we used the hydrolysis of Ub- 
AMC here to test the deubiquitinating activity in vitro. 
The kinetic curves of proteases and variants with the 
substrate Ub-AMC were hyperbolic and the initial velo- 
cities were not linear over the concentration of substrate. 
However, saturation was still not observed within a reas- 
onable time. Only when the ratio of protein to substrate 
was 1:1 or larger, were we able to achieve saturation 
within a limited time (data not shown). As the initial ve- 
locities did not increase in a strictly linear fashion with 
substrate concentration, application of the equation 


WIE] tot, = KapplS] to mimic k,,,/Ky, would lead to large 


standard errors. We were however able to fit the data to 
the Michaelis-Menten equation using the GraphPad 
Prism program (GraphPad Software), even though satur- 
ation could not be reached (Supplemental Figure 1). All 
assays were performed in duplicates. 


RESULTS 


Overall structure of MERS-CoV PL”” in complex 
with human ubiquitin 

The substrate-binding site of MERS-CoV PL” features 
significant differences from those of the corresponding 
SARS-CoV enzyme and human ubiquitin-specific pro- 
teases (USPs, such as, USP14) (Hu et al., 2005; Chou et 
al., 2014; Ratia et al., 2014). It is therefore of interest to 
determine the crystal structure of the complex between 
MERS-CoV PL”” and its substrate, human ubiquitin. 
Hence, we crystallized the ubiquitin (Ub) complex of a 
MERS-CoV PL”” variant that had the active-site Cys111 
replaced by serine (C111S) and determined the structure 
at 3.16 A (Figure 1A). 

There is one PL” (C111S)-Ub complex per asymmet- 
ric unit. Using the PDBePISA server (Krissinel & Hen- 
rick, 2007), the total interface region of MERS-CoV 
PL”® (C111S)-Ub was determined as 813 A’, close to the 
915 A’ interface of SARS-CoV PL”’-Ubal (ubiquitin al- 
dehyde; PDB entry: 4MM3, Ratia et al., 2014) and the 
999 A* of SARS-CoV PL” (C112S)-Ub (PDB: 4M0W, 
Chou et al., 2014), but less than the 1503 A’ observed for 
USP14-Ubal (PDB: 2AYO, Hu et al., 2005). The overall 
structure of MERS-CoV PL” (C111S) is very similar to 
that of the substrate-free PL”® (Lei et al., 2014), with a 
root-mean-square deviation (RMSD) of 0.91 A for cor- 
responding Ca atoms between these two structures. Nev- 
ertheless, two differences are immediately visible: 1), the 
zinc-finger motif has moved and closed in onto the Ub, 
compared to the free PL”’. The zinc ion position has 
shifted by about 4 A and the largest deviation between 
the two structures is ~6 A for Cys228 of the zinc-finger 
region (Figure 1B); 2) the mobile loop *"GIETAVG 
also named “BL2 loop”, is defined by clear main-chain 
electron density (Figure 1C). This loop is disordered in 
substrate-free PL” (Lei et al., 2014; Bailey-Elkin et al., 
2014). Compared to MERS-CoV PL?”, the BL2 loop 
*’GNYQCG” is shorter by one residue in SARS-CoV 
PL™. 

At the time when we deposited in the PDB the co- 
ordinates for the crystal structure of the non-covalent 
complex MERS-CoV (C111S)-Ub (PDB: 4WUR), 
Bailey-Elkin et al. (2014) described two crystal struc- 
tures for a covalent complex MERS-CoV PL””°-Ub, in 
which an alkyl bromide group introduced at the C-ter- 
minus of Ub had formed a thioether with the active-site 
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Ub! domain 


Figure 1. Structure of MERS-CoV PL”” in complex with 
human ubiquitin (Ub). (A) Cartoon view of the overall 
complex structure. PL’® is shown in green and Ub in 
red. The N and C termini are marked by underlined let- 
ters, and the PL’ subdomains are divided by black 
dashed lines. The ubiquitin-like (Ubl) domain, thumb, 
palm, and fingers subdomains are indicated. The cata- 
lytic triad Cys111-His278-Asp293 is shown by spheres 
(yellow, blue, and red; Cys was replaced by Ser in this 
study). The zinc atom is displayed as a gray sphere. 
Five PL” regions interacting with Ub are colored dark 
blue and labeled with Roman numbers (I-V). (B) Super- 
position of the fingers domain of Ub-bound MERS-CoV 
PL°® and the substrate-free enzyme (PDB: 4P16; Lei et 
al., 2014). The zinc atom is shifted by ~4 A, and the 
largest difference is ~6 A. (C) The BL2 loop (blue; *”'GI- 
ETAVG’”’) of PL”. A 2F.-F, electron density (gray; 1.0 
0) is displayed. The side-chains of residues in this loop 
are only partly defined by electron density. 


Cysl11 of MERS-CoV PL?”. These structures were in 
space groups P6,22 and P6, and were named “closed” 
and “open” PL”’-Ub complexes (PDB: 4RFO and 4RF 1), 
respectively (Bailey-Elkin et al., 2014). Even though the 
zinc-finger motif of PL” in the former complex (space 
group P6,;22) is closer to the Ub than in the P6, complex, 
both the “closed” and “open” complex show almost the 
same interactions between PL” and Ub (Bailey-Elkin et 
al., 2014). The differences may be caused by crystal 
packing. Our non-covalent complex is similar to the 
“open” form of the covalent PL”’-Ub complex (RMSD = 
0.58 A for the PL”” and 0.85 A for Ub, based on all Ca 
atoms). 

All parts of the PL” except for the ubiquitin-like 
(Ub1l) domain interact with Ub. Most interactions in- 
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Figure 2. Interactions between MERS-CoV PL®”® and 
Ub. (A) Five regions of PL"® bind to two areas of Ub. 
The surface of PL” is shown in wheat color. The five 
interacting regions are shown in brown, cyan, yellow, 
green, and purple, and they are also labeled by Roman 
numbers. The two interacting areas of Ub are marked 
“A” and “B”. Region A is highlighted by dark blue dots, 
and region B (the RLRGG sequence; Ub residues are in 
italics) is depicted in the ball-and-stick style. (B) Region 
A of Ub interacting with PL’®. A cartoon view of region 
A is shown in purple, and $3 and £4 are labeled. 
Residues of Ub (purple) and PL” (green) are dis- 
played in the ball-and-stick style, and labeled in red and 
black, respectively. The 2F,-F, electron density (blue; 
1.0 0) of the side-chains of Arg42 and Asp165 is dis- 
played. Hydrogen bonds are indicated by black dashed 
lines, and the salt-bridge between R42 and D165 is de- 
picted by two red dashed lines. (C) The RLRGG bind- 
ing site of PL"°. The P5—P1 residues are shown in the 
ball-and-stick style. Regions |, II, IV, and V that interact 
with the P5—P1 residues are labeled. 


volve five surface regions of the PL” and two regions of 
Ub (Figure 2A). These five regions of the PL”” are la- 
belled by Roman numbers: I, Leul06—Tyr1 12; II, 
Alal62-Arg168; III, Cys208—-Val210; IV, 
Gly248—Pro250; V, Phe269-Tyr279. I and IJ are situ- 
ated in the thumb domain; III is in the fingers domain; 
and IV and V are in the palm domain (Figures 1A and 
2A). In the following and in the figure labels, residues of 
Ub are indicated in italics to distinguish them from 
residues of PL””. The two interacting regions of Ub are: 
A, Arg42—Gin49; B, Arg72-Gly76 (Figure 2A). Region 
A of Ub consists of 83, 64, and the loop between B3 and 
B4 (Figure 2B). Regions II and III of PL”” interact with 
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region A of Ub. Region B comprises the five C-terminal 
residues, RLRGG, and is in contact with regions I, II, IV, 
and V of PL” (Figure 2C). The C-terminal RLRGG mo- 
tif contributes the majority of the interactions with the 
PL”; the buried surface between Ub region B and PL”” 
is 477 A* (out of a total of 813 A’). In order to make 
these interactions clear, we describe here in some detail 
the contacts between the PL”” and regions A and B of Ub. 


Interactions of MERS-CoV PL”” (C111S) with Ub 
region A 

Region A (Arg42—Gln49) of Ub inserts into the space 
between the thumb and fingers domains of PL”” (Fig- 
ures 2A—B). Residues [/e44, Ala46, and Gly47 engage in 
hydrophobic interactions with Tyr209 and Val210 of re- 
gion III of PL” (Figure 2B). The hydrophobic patch 
(1le44, Ala46, and Gly47) is a common interaction re- 
gion utilized by Ub-binding proteins (Dikic et al., 2009). 
In particular, the interaction of /e44 with Val210 of 
PL” is important for the deubiquitinase (DUB) but not 
for the protease activity. The variant V210R shows dra- 
matically reduced DUB activity, as demonstrated by 
Bailey-Elkin et al. (2014) (according to the numbering 
scheme of these authors, V210 is V1691). 

A salt-bridge exists between region A of Ub and the 
PL” (Figure 2B), namely between the side-chains of 
Arg42 and Asp165 (region II in PL”’). In addition, a hy- 
drogen bond is formed between the main-chain O atom 
of Gly47 and the main-chain amide of Val210 (region 
III). In the SARS-CoV PL™°(C112S)-Ub complex (Chou 
et al., 2014), Arg42 forms a salt-bridge with the negat- 
ively charged Glu168, a residue which is replaced by the 
positively charged Arg168 in MERS-CoV PL””. Con- 
sequently, the same salt-bridge cannot be formed in the 


MERS-CoV PL?” (C111S) -Ub complex. Instead, 
MERS-CoV PL”” has Asp165 interacting with Arg42. In 
our deubiquitinating (DUB) kinetic assays, the DI65A 
variant shows a dramatically reduced DUB activity; the 
K.,/Ky is about 78-fold decreased compared to that of 
the wild-type (Table 2). The Km value of the D165A 
variant is about 4-fold higher than that of the wild-type 
MERS-CoV PL?”, suggesting that this amino-acid re- 
placement reduces the Ub binding affinity. Meanwhile, 
the D165E amino-acid replacement shows a catalytic ef- 
ficiency comparable to wild-type towards Ub-7-amino-4- 
methylcoumarin (Ub-AMC) (Table 2). However, we no- 
ticed that the K,, value of D165E for the DUB activity is 
about 2-fold larger than for the WT enzyme. Although 
Glul65 can mimic Asp165 here, we propose that the 
longer side-chain of Glu may fit less perfectly compared 
to Asp. These results indicate that the salt-bridge 
between Asp165 and Arg42 could be important for the 
PL’’’’s DUB activity, in addition to the interaction 
between Asp165 and the P4-amide group (see below). 


Interactions of MERS-CoV PL””® (C111S) with Ub 
region B 

Region B of Ub comprises the five C-terminal residues, 
RLRGG. These five residues bind to the narrow active- 
site channel between the thumb and palm domains of the 
PL” (Figure 2C). They mainly interact with regions I, II, 
IV, and V of the protease. Residues RERGG are compat- 
ible with the PL”” cleavage motif, (R/K) (L/I) XGG 
(P5—P1), in the MERS-CoV polyproteins; their interac- 
tions with the PL”” are discussed here in terms of sub- 
sites S1 to SS. 


SI and S2 subsites. Pro163 (region II of PL”) and the 


Table 2. Kinetic parameters of MERS-CoV PL®® and SARS-CoV PL°® 


Z-RLRGG-AMC Z-LRGG-AMC Ub-AMC 

Kapp (umol min”')* Ky (umol/L) Keat (min’') KeatlKy (umol/L“ min”) 
MERS-CovV PL” 
WT*' (1.3+0.1)x10° (1.040.01)x10° 6140.7 8.4+0.7 1.4+0.2 
F269Y (2.0+0.1)x10° (1.440.2)x10° 85+0.1 22.44 1.5 26202 
D165E (ee 2) 108 (Gad) x0n lo eee 20 £102 7203 
D165A (8.2+1.6)x10° (1.0+0.2)x 10% 23.9 + 6.2 0.42 + 0.04 0.018 + 0.005 
D164E (1.540.1)x10° (8641.8)x10* 15.7444 (Eicies oe 1.0+0.3 
D164A (2.940.1)x10* (6.5+0.8)x10* 6540.6 2.5 +0.01 0.40 + 0.04 
wi"? (3.0 + 0.1) x 10° 1432 20 18.8 + 1.2 (ee 
SARS-CoV PL’” 
wT! 0.23 + 0.04 (3.540.5)x10°% 57.2+10.9 197.4 + 28.0 G50 
wi"? oof 01 50.6+7.4 75.9 + 8.1 |Se02 


Notes: *: Wild type; * : Kapp 


is an approximation to k,o/Kjy; ': Our results (this work); *: Baez-Santos et al., 2014b. 
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side-chains of Asn109 and Tyr112 (located in region I) 
form a space-restricted S1 site to accommodate Gly76 
(P1). The carbonyl oxygen atom of Gly277 (region V) 
accepts a hydrogen bond from the amide of Gly76 (Fig- 
ure 3). 

The side-chains of Tyr112 (region I) and Phe269, 
Val276 as well as Tyr279 (region V), and Gly277 form a 
restricted space for Gly75 (P2). Two hydrogen bonds 
link the main-chain at Asp164 of the PL”” and Gly75 
(Figure 3). The S1 and S2 sites of the protease are too 
small to accommodate any other residue but glycine. 


S3 subsite. The main-chain O atom of Arg74 (P3) ac- 
cepts a hydrogen bond from the Gly277 amide. In the 
complexes SARS-CoV PL”*(C112S)-Ub or SARS-CoV 
PL”’-Ubal, the main chain at Arg74 forms two hydrogen 
bonds (Chou et al., 2014; Ratia et al., 2014), namely with 
the amide of Gly272 and with the hydroxyl group of 
Tyr265. The former hydrogen bond is conserved in the 
MERS-CoV PL”” (C111S)-Ub complex, but the latter is 
not. Tyr265 of SARS-CoV is replaced by Phe269 in 
MERS-CoV, which lacks the ability to form a hydrogen 
bond with the main-chain amide of Arg74. In agreement 
with this difference, the DUB activity of SARS-CoV 


R168 


Figure 3. Details of the interactions between the C-ter- 
minal RLRGG of Ub and MERS-CoV PL®”. RLRGG 
residues are shown in purple in the ball-and-stick style, 
and they are labeled in red. For clarity, the 2F,-F, elec- 
tron density (blue; 1.0 o) of the RLRGG main chain is 
shown. PL” residues are displayed in green in the ball- 
and-stick style, and they are labeled in black. Hydrogen 
bonds are displayed as black dashed lines, and the 
salt-bridge between R72 and D164 is depicted as two 
red dashed lines. 
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PL” is 2.5-fold higher than that of MERS-CoV PL” in 
our kinetic assay; furthermore, the MERS-CoV PL” 
F269Y amino-acid replacement leads to enhancements 
by about 1.5-, 1.4-, and 1.9-fold of the hydrolytic activit- 
ies towards carbobenzoxy-Arg-Leu-Arg-Gly-Gly-7- 
amino-4-methylcoumarin (Z-RLRGG-AMC), Z-LRGG- 
AMC, and Ub-AMC, respectively (Table 2). 

The side-chain of Arg74 is exposed to the solvent in 
our MERS-CoV PL’” (C111S)-Ub complex. This side- 
chain shows remarkable variability in the interactions it 
makes in the different complexes. In the “open” but not 
in the “closed” covalent MERS-CoV PL”’’-Ub complex 
(Bailey-Elkin et al., 2014), it donates a hydrogen bond to 
the main-chain carbonyl oxygen of Thr1755 of the BL2 
loop (corresponding to Thr274 in our numbering 
scheme). In SARS-CoV PL”°-Ubal (Ratia et al., 2014) 
but not in the SARS-CoV PL”” (C112S)-Ub complex 
(Chou et al., 2014), the side-chain of Arg74 forms a hy- 
drogen bond with the main-chain carbonyl oxygen of 
Gln270. Instead, Arg74 is involved in a relatively weak 
salt-bridge with Glul62 in the SARS-CoV PL’”® 
(C112S)-Ub complex (Chou et al., 2014). None of these 
interactions exist in our MERS-CoV PL”’-Ub complex. 


S4 subsite. The main-chain amide of Leu73 (P4) donates 
a hydrogen bond to the side-chain of Asp165 (region IT) 
(Figure 3). This hydrogen bond is conserved in all the 
complexes compared here. In our peptide-cleavage as- 
say, the D165A variant shows about 16-fold and 10-fold 
lower activities towards substrates Z-RLRGG-AMC and 
Z-LRGG-AMC, respectively, indicating that Asp165 is 
not only important for interacting with Arg42 of Ub (see 
above). 

The side-chain of Lew73 is embedded in a hydrophobic 
pocket which is formed by the CB atom of Asp165, the 
side-chain of Pro250 (region IV), Phe269, as well as the 
CB and Cy atoms of Glu273 (region V). Asp165 and 
Pro250 are conserved in SARS-CoV PL”” (Asp165 and 
Pro249). Phe269 is replaced by Tyr265, and Glu273, 
situated in the BL2 loop, is replaced by Tyr269. 
However, the side-chains of all these non-conserved 
residues possess the ability to provide a hydrophobic en- 
vironment to accommodate Leu73. 


S5 subsite. The side-chain of Arg72 (P5) is located 
between the PL” thumb domain and region A of ubi- 
quitin. It forms a salt-bridge with the side-chain of 
Asp164 in MERS-CoV PL?”. In addition, the guanidin- 
um group of Arg72 may be involved in a 2-7 interaction 
with that of Arg168 (Figure 3). These interactions have 
also been described for the covalent complex of MERS- 
CoV PL” with Ub (Bailey-Elkin et al., 2014). Arg72 is 
not subject to strict space limitations; accordingly, this 
residue displays different binding patterns in the two 
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SARS-CoV PL”’-Ub complexes. In the covalent SARS- 
CoV PL”’’-Ubal complex (Ratia et al., 2014), Arg72 
forms a salt-bridge with Glu168. In the non-covalent 
SARS-CoV PL” (C112S)-Ub complex (Chou et al., 
2014), Arg72 is exposed to the solvent and does not in- 
teract with Glu168 (Arg42 instead forms a salt-bridge 
with Glul68, as mentioned above). In our kinetic assay, 
the D164A variant of MERS-CoV PL” displays a ~4.5- 
fold and a ~3.5-fold reduced activity, respectively, for Z- 
RLRGG-AMC and Ub-AMC (Table 2). For Z-LRGG- 
AMC, the activity is decreased just a little (by about 
~1.5-fold), because there is no P5-Arg in this substrate 
(Table 2). These data demonstrate that Asp164 is import- 
ant for the interaction with Arg72. 

In summary, the main-chain heteroatoms of P5—P1 
form a hydrogen-bonding network with PL”. The bind- 
ing characteristics of P1, P2, and P4 are conserved in all 
MERS-CoV and SARS-CoV PL”’’-Ub complexes. 
However, P3-Arg and P5-Arg assume binding patterns 
that differ between the various MERS-CoV and SARS- 
CoV PL”’-Ub complexes. 


DISCUSSION 


Viral proteins are likely to possess multiple functions, as 
exemplified by non-structural protein 1 (NS1) of influ- 
enza A viruses (Hale et al., 2008), the nucleocapsid (N) 
protein of coronaviruses (Chang et al., 2014), or the 
Nsp14 exonuclease-guanyl-7-methyltransferase of 
coronaviruses (Chen et al., 2009). Exhibiting DUB and 
proteolytic activities, the papain-like proteases of 
coronaviruses are no exception here. Although the over- 
all folds are conserved, the enzyme activity and sub- 
strate-binding modes of CoV PL”’’s differ in detail. 
Therefore, no coronavirus PL”” can be considered a gen- 
eral model for all its homologues (Baez-Santos et al., 
2014b). 


MERS-CoV PL”” is not optimized for catalytic 
efficiency 

We have previously noticed that the oxyanion hole of the 
MERS-CoV PL?” appears to be deficient (Lei et al., 
2014; also see the discussion below). Similarly, the re- 
cognition of ubiquitin by the enzyme appears to be sub- 
optimal. Thus, the main-chain amide of the P3-Arg 
residue has no hydrogen-bonding partner on the MERS- 
CoV PL””, because the near-by side-chain of Phe269 is 
incapable of accepting an H-bond. The corresponding 
residue is Tyr265 in the SARS-CoV PL”, which is per- 
fectly positioned to accept the hydrogen bond from the 
P3-amide. Indeed, when we replaced Phe269 by Tyr in 
MERS-CoV PL?”, the DUB activity of the enzyme in- 
creased by a factor of almost 2 and the peptidolytic activ- 


ity by ~1.5. The evolution of viral enzymes is not neces- 
sarily driven by optimization of catalytic efficiency. This 
is particularly true for viral proteases that have to ensure 
the availability of non-structural proteins in the correct 
temporal order when they cleave them out of the viral 
polyproteins; in fact, too rapid a polyprotein processing 
might be counterproductive. On the other hand, a more 
efficient DUB activity should help the virus in counter- 
acting the innate immune response of the host cell. As 
the binding of the P5—P1 residues of ubiquitin and of the 
polyprotein cleavage sites obviously influences both the 
DUB and proteolytic activities of the PL””, the subop- 
timal catalytic efficiency that we observe may be a con- 
sequence of a compromise between the requirements of 
the two activities. 

With regard to the oxyanion hole, we previously pro- 
posed that the backbone amide of Asn109 (located in a 
B-turn connecting B7 and a4, Figure 4) may contribute to 
the stabilization of the oxyanion intermediate in PL”” 
catalysis, along with the main-chain amide of the active 
site Cys111, although this may require a slight rearrange- 
ment of this B-turn (Lei et al., 2014). In our non-cova- 
lent MERS-CoV PL””-Ub complex, we do not see any 
rearrangement of the B-turn. More or less in agreement 


Figure 4. The side-chain of Asn109 is unlikely to con- 
tribute to the oxyanion hole of MERS-CoV PL”®. PL°” is 
shown in green in cartoon view. N and C termini are 
marked by underlined letters. Strand B7, helix a4, and 
the loop between them are colored in blue. Asn109 and 
Gly161 are indicated in ball-and-stick style in cyan and 
purple, respectively, and they are labeled in black. The 
hydrogen bond between them is displayed as a black 
dashed line. The active-site Cys111 (replaced by Ser in 
this study) is shown in green in ball-and-stick style. 
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with our suggestion, Bailey-Elkin et al. (2014) proposed 
on the basis of their covalent MERS-CoV PL”’-Ub com- 
plex structure that the main-chain amides of Asn1590, 
Asn1591, and Cys1592 (corresponding to Asn109, 
Asn110, and Cys111 in our numbering scheme) form the 
oxyanion hole. On the other hand, Lee et al. (2015) ar- 
gued that the side-chain of Asn109 could contribute to 
the oxyanion hole. They found that the N109A replace- 
ment completely abolished enzyme activity. As we re- 
ported earlier (Lei et al., 2014), the side-chain amide of 
Asn109 makes a strong hydrogen bond with the con- 
served Glyl61 (Figure 4). Any reorientation of the 
Asn109 side-chain towards the oxyanion would require a 
disruption of this strong interaction; this is not very 
likely. In conclusion, lacking a side-chain in the proper 
spatial orientation and capable of donating a hydrogen 
bond to the oxyanion (such as Trp107 in SARS-CoV 
PL”), the oxyanion hole of MERS-CoV PL”” seems to 
be deficient. 


Unique features of Ub recognition by MERS- 
CoV PL’? 

Apart from the less than optimum binding of the P3 
residue, there are other differences in the way MERS- 
CoV PL’”® and SARS-CoV PL”” recognize human ubi- 
quitin. The formation of a salt-bridge between Arg72, the 
P5 side-chain of Ub, and Asp164 is unique for MERS- 
CoV PL”. The same interaction exists in the two cova- 
lent MERS-CoV PL”’-Ub complexes (Bailey-Elkin et 
al., 2014). This binding mode is very different from the 
Glul168 - Arg72 salt-bridge in the SARS-CoV PL””-Ubal 
complex (Ratia et al., 2014). Glu168 is conserved in 
HCoV-NL63 PL2”” and replaced by Asp in the HCoV- 
229E, HCoV-OC43, and HCoV-HKUI1 PL2”"s (for se- 
quence alignments, see Barretto et al., 2005; Baez-San- 
tos et al., 2014b), so the same type of interaction is likely 
to be realized in the Ub complexes of these enzymes. 
However, the corresponding residue in MERS-CoV PL”” 
is Arg168; hence, Asp164 is used instead to bind Arg72. 
Asp 164 is in fact unique in MERS-CoV. It is replaced by 
Gly in SARS-CoV PL” and the PL2”"s of HCoV NL63 
and HCoV 229K, by Ala in HCoV-OC43 PL2””, and by 
Ser in HCoV-HKU1 PL2””. Our kinetic results for the 
D164A replacement (see Results) emphasize the import- 
ance of the unique Asp164 residue. 


The role(s) of the ubiquitin-like (Ubl) domain 

In our PL” (C111S)-Ub complex, the Ub! domain of 
PL” shows no interaction with ubiquitin. The relative 
orientation of the Ub! domain in the substrate-bound 
PL” is the same as in substrate-free PL””. The Ubl do- 
main of MERS-CoV PL”” is not required for the IFN 
antagonism activities (Baez-Santos et al., 2014b), but it 
is required in case of SARS-CoV PL”” (Frieman et al., 
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2009). Recently, Mielech et al. (2015) reported that the 
Ub! domain of mouse hepatitis virus (MHV) is an im- 
portant modulator of PL” stability and viral pathogenes- 
is. Although the Ubl domain shows variable effects in 
different CoVs, the high degree of conservation of the 
domain throughout the family suggests that it may play a 
common biological role. One possible function is that the 
Ub! might be involved in protein-protein interactions. 
Ubiquitin-like domains are known to function as binding 
modules in such interactions. For example, the kinase 
Raf contains a Ub! domain for interaction with human 
Ras (Fetics et al., 2015). Also, human ubiquitin-specific 
protease (USP) 7 includes five Ubl (1-5) domains, of 
which the second is bound by the Herpes simplex virus-1 
immediate-early protein ICPO to counteract the intrinsic 
antiviral response of the host cell (Pfoh et al., 2015). 


PL”” inhibitors 

The BL2 loop of SARS-CoV PL”” is important for bind- 
ing inhibitors (Baez-Santos et al., 2014a; Lee et al., 
2015). This loop is variable in different CoV PL’’s. A 
potent inhibitor of the SARS-CoV PL””, N-[(3-fluoro- 
phenyl) methyl]-1-[(1R)-1-naphthalen-1-ylethyl] piperid- 
ine-4-carboxamide (compound 3k, IC.) = 0.15 + 0.01 
umol/L) was found to have no effect on the MERS-CoV 
enzyme (Baez-Santos et al., 2014b). Tyr269 and GIn270 
of SARS-CoV PL””, which are important for binding this 
inhibitor (Baez-Santos et al., 2014a), are replaced by 
Glu273 and Ala275 in the MERS-CoV protease. It seems 
that this structural difference in the BL2 loop has a re- 
markable impact on the effectiveness of the inhibitor. 
The structure of the MERS-CoV PL” (C111S)-Ub com- 
plex presented here will facilitate virtual screening of 
chemical libraries for specific anti - MERS-CoV PL” 
inhibitors (Hilgenfeld, 2014). 


The di-Ub site of PL’ 

SARS-CoV and MERS-CoV PL?’s can digest K48- and 
K63-linked (poly)ubiquitin chains and remove ISG15 
from proteins covalently linked to it (Ratia et al., 2014, 
Baez-Santos et al., 2014b). SARS-CoV PL”” prefers 
binding of K48-Ub, and ISG15 over mono-ubiquitin 
(Ratia et al., 2014). All evidence suggests that on the 
PL””, at least two ubiquitin-binding sites (or a binding 
site for diubiquitin-like molecules such as ISG15) exist. 
Ratia et al. (2014) proposed hypothetic models for com- 
plexes of SARS-CoV PL” with K48-Ub, and ISG15. 
These authors identified two major hydrophobic binding 
sites on the PL”” for the first Ub (Ub1) and the second 
Ub (Ub2). The binding site for Ub] comprises Met209, 
Pro248, and Pro249. This hydrophobic patch is con- 
served in MERS-CoV PL””, although the residues are 
not exactly the same. Met209 of SARS-CoV PL” is re- 
placed by Val210 (region IIIT) in MERS-CoV. Pro248 is 
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replaced by Thr249, whereas Pro249 is conserved 
(Pro250; region V) (Figure 5). The hypothetic second 
binding site for Ub2 (also named “ridge” region, Ratia et 
al., 2014) is located to the first o helix (a2) in the thumb 
domain, including residues Phe70, His74, and Leu76 of 
SARS-CoV PL”. According to a structural alignment of 
the SARS-CoV and MERS-CoV PL”’’s, Phe70, His74, 
and Leu76 are changed to Lys69, Gly73, and Val75, re- 
spectively (Figure 5). In SARS-CoV PL””, the F70S and 
F70A replacements lost the affinity to K48-Ub, and 
ISG1I5 in vitro (Ratia et al., 2014); therefore, the pres- 
ence of Lys69 in MERS-CoV PL”” instead of Phe70 
strongly suggests that this enzyme should bind Ub2 ina 
different way. Bailey-Elkin et al. (2014) predicted 
Asn1673 and Vall1674 of the fingers subdomain (corres- 
ponding to Asn192 and Val193 in our PL”) as the distal 
Ub site of K63 di-Ub, but they found that their DUB 
activity data do not support their prediction for the Ub2 
binding site. The structure of the MERS-CoV PL’”-Ub 
complex reported here reveals the exact Ub! binding site 
on the PL”, but the Ub2 binding site should be identi- 
fied by crystallizing the PL” in complex with poly- or 
di-ubiquitin. 

In summary, the crystal structure of the MERS-CoV 
PL”’-Ub complex provides valuable information that 
helps understand the multiple functions of coronavirus 
papain-like proteases. Mutational studies additionally 
highlight features of the MERS-CoV PL””. The different 
substrate-binding patterns should be kept in mind when 
designing inhibitors for PL”’s of CoVs, even though the 
overall structures of these enzymes are conserved. Fur- 
thermore, it would be helpful to obtain the structure of 


‘proposed 
Ub? site’, 


Figure 5. The two Ub-binding sites of PL". The Ub1 
and the proposed Ub2 (according to a structural align- 
ment with the SARS-CoV PL®®: Ratia et al., 2014) bind- 
ing sites are depicted as purple dots. The N and C ter- 
mini of PL®’° are marked by underlined letters. All 
residues related to the two Ub binding sites are labeled. 
The catalytic triad Cys111-His278-Asp293 is indicated 
by yellow, blue, and red spheres. 
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PL” in complex with di-Ub or ISG15 in the future. 
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